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Abstract 

The subprocess gg — ► tc + ic in the two-Higgs-doublet 
model with flavor- changing scalar couplings is exam- 
ined at the one loop level. With perturbative QCD 
factorization theorem, the corresponding cross sections 
for hadron-hadron collisions are computed numerically. 
The results are applicable to the whole mass range of 
the weakly coupled Higgs bosons. In case we could effi- 
ciently exclude the severe backgrounds of the tc{tc) pro- 
duction signal, probing the flavor- changing top-charm- 
scalar vertex at hadron colliders would be very promis- 
ing and accessible experimentally 
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I. Introduction 

There are stringent experimental constraints against the existence of tree level neutral flavor 
changing interactions, especially for light quarks [0]. This leads to the naturalness to suppress the 
flavor changing neutral current interactions (FCNCs) for all kinds of model building studies, that is 
realized in terms of the Glashow-Iliopoulos-Maiani (GIM) mechanism in the standard model(SM). 

The two-Higgs-doublet models(THDMs) are extensions of the SM in which one more scalar dou- 
blet is added. In order to forbid possible tree level FCNCs to appear in the models, Glashow and 
Weinberg proposed a neutral flavor conservation (NFC) condition by imposing discrete symmetries 
on the models. The common THDMs with the NFC condition can be divided into two categories, 
i.e., Model I and Model II. In Model I, both the up- and down- type quarks couple to the same one 
of the Higgs doublets, but in the Model II the up- and down-type quarks couple to the two Higgs 
doublets respectively. In fact, the suppression facts of the FCNC processes in the cases involving 
the light down-type quark sector have been observed experiment ally [|], whereas those relevant to 
the up-type quark sector have not been established so well. Therefore the so-called Model-Ill of 
the THDM is proposed and it allows the FCNC Yukawa couplings which has the character that the 
couplings are related to the masses of the coupled flavors even at tree level0. 

As pointed out by Cheng, Sher and other authors 00 H@, the Yukawa couplings are typically 
proportional to the masses of the coupled fermions at the vertices, it is rather natural to imagine 
having such Yukawa couplings for the FCNC interactions instead of placing the constraints on the 
theory. In this case, low energy limits on the FCNCs may be evaded, because the flavor changing 
couplings to the light quarks are small, and the suppression of the FCNC processes involving light 
quarks can be automatically satisfied. Then the imposition of discrete symmetries to obtain the 
NFC condition, therefore, is unnecessary, which is normally invoked in common THDMs to prevent 
the FCNC interactions at tree level. In the literature such the THDM without the NFC condition 
is called as the 'third' THDM (i.e., THDM III) |§, where the up-type and down-type quarks are 
allowed simultaneously to couple to more than one scalar doublet. In the framework of the model 
the effects of the FCNC interactions involving the heavy quark will be enhanced. 

Now the top quark has been observed and a unexpected very large mass m t = 173.8 ± 5.2 GeV 
(world average value) has been recognized 0. This 'extraordinary' mass scale of the top quark may 
have many important implications pertaining to many outstanding issues. One of them is to test 
the FCNC processes requested in the Model-Ill. Namely, it is a 'good place' to observe the facts 
of the Model-Ill for the FCNC interactions. The measurement of FCNC processes involving top 
quark would provide an important test for the discrimination of various models. In the THDM 
III one would expect that large effects of the FCNC interactions could manifest themselves in the 
cases involving the massive top quark. Therefore testing the existence of the flavor changing scalar 
interactions involving top quarks are clearly important. D. Atwood et al. [JTO]] presented results 
of a calculation for the process e+e~ -> tc(or tc)in the THDM III, and they obtained i? tc /A 4 to 
be in the order of 10~ 5 with proper parameters. Recently, Jiang et al.|TT[ studied the production 
rates of e + e~ — > 77 — > tc + ic for the NLC and found that this process is more promising than the 
straight production via e + e~ collisions for probing the FCNC interactions. Abraham et al.[T^] also 
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investigated the anomalous tqj couplings via single top quark production processes by considering 
effective Lagrangian of the lowest dimension with 77 collisions, and found the processes would be 
observable with suitable strength of anomalous coupling as long as b-tagging and suitable kinematic 
cuts are taken properly. 

After the termination of the running of the LEP2, the hadron colliders Tevatron and LHC will 
be the only machines in searching for the FCNC processes. It is believed that more experimental 
events involving top quark will be collected in these hadronic machines. It will give a good chance 
to study the physics relevant to the FCNC processes of the Model-Ill. In this paper we are to 
study the problem and present complete one-loop calculation for the subprocess gg — > tc(or tc) to 
the order 0[m t m c j m^) for the THDM Model-Ill. In fact, the obtained results in the paper are 
applicable to the whole mass range for weakly coupled Higgs bosons. The production cross sections 
of pp(p) — > gg —>■ tc + tc + X are also given for the Tevatron and LHC. The paper is organized as 
follows: The details of the calculation are given in Sec. II. In Sec. Ill numerical results, discussions 
and a short summary are presented. Finally, the explicit expressions used in the paper are collected 
in appendix. 

II. Calculation 

In the third type of the two-Higgs-doublet model, the up-type and down-type quarks are allowed 
simultaneously to couple to more than one scalar doublet. Since there is no global symmetry that 
distinguishes the two doublets in the model, we will assume that only one of them ((pi) develops a 
vacuum expectation value and the second one (^2) remains unbroken, i.e. 

< 0i >= ( ° u/v /2 ) , <0 2 >=O (1) 

where v ~ 246 GeV. The physical spectrum of Higgs bosons consists of two scalar neutral bosons h° 
and H°, one pseudoscalar neutral boson A and two charged Higgs H ± , 

H° = ^/2[(Re<j)\ - v) cos a + sin a] , 

h° = y/2[- (Re4>l - v) sin a + Re}° 2 cos a] , (2) 
A = s/2(-Im(j>l). 

The masses of the five neutral and charged Higgs bosons and the mixing angle a are free parameters 
of the model. The Yukawa couplings to quarks are [[J, 

C% = XfjQifaUj + \',]Q ,<:>: Pj + gjQifaUj + CllQ.oJ),. (3) 

where the first two terms give masses of the quark mass eigenstates, and £^ and £E are the 3x3 
matrices which give the strength of the flavor changing neutral scalar vertices. The £s are all free 
parameters and can be constrained by the experimental data. If we neglect CP violation, the £s are 
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all real. We will use the Cheng-Sher Ansatz(CSA)[K| and let 



«i ~ " 



And we can parametrize the Yukawa couplings as 



= g— — — ^-A. (4) 

mw 

Comparing it with the usual gauge couplings of SU{2) x U(l), one has A = In our calculation 
we use A = -j= and note that there is no theoretical bound on the coupling factor A. 
The subprocess producing tc(tc) via gluon-gluon collisions, 

gg -> tcitc) 

can be induced through one-loop diagrams at the lowest order, and the Feynman diagrams are drawn 
in figure 1(a) and figure 1(b), where the contributions from the one-loop diagrams involving neutral 
and charged Higgs bosons are presented, respectively. The contributions from the figures of Fig. 1(b) 
with the charged Higgs boson in loops being replaced by W-boson, is much smaller due to the GIM 
suppression and Yukawa coupling. We can omitted this part in our calculation. The diagrams 
exchanging the two external gluon-gluon lines are not shown, but are numbered in Fig. 1(a) and 
Fig. 1(b). Fig.l(a)(l ~ 12) and Fig. l(b)(l ~ 6) are the self-energy diagrams, Fig. l(a)(13 ~ 20) and 
Fig. l(b)(7 ~ 10) are the vertex correction diagrams, Fig. 1(a) (25 ~ 34) and Fig. 1(b) (13 ~ 15) 
are the s-channel diagrams, Fig. l(a)(21 ~ 24) and Fig. l(b)(ll ~ 12) are the box diagrams. Note 
that in the present case at one-loop level the ultraviolet divergence would be canceled automatically, 
if all the one- loop diagrams at the 0{m t m c I my) order in the THDM III are included. In this work, 
we perform the calculation in the t'Hooft-Feynman gauge. 

To simplify the calculation we set a = and let all scalar bosons be degenerate, i.e., m^o = m^o = 
rriH± = M s where M s is the common scalar mass. The contribution from the coupling involving H° 
is suppressed due to a = 0. 

In the calculation for the s-channel diagrams(Fig.l.(a)(25 ~ 28)), we take into account the width 
effects of the h° and A propagators. The decays of h° to WW and ZZ are suppressed, because of 
the factor sin a in the h°WW and h°ZZ couplings, and h° decay to A A is also forbidden due to 
the case of the degenerate masses of h° and A . Note that the pseudoscalar A does not couple with 
gauge boson pair. Therefore only the decays of h° and A to final states need to be considered, 
where qi and qj represent quarks of flavor i and j respectively. The decay width for the scalar h° can 
be written as[[L4]] 



r(/^ 9 ?) = ^ m »(i-4 m ?«„) 3/2 



and 



3g 2 m h a 



T(h u ^tc + tc) = „ • 2m t m c l-{m t + m c ) /m h0 x l - (m t - m c ) /m h0 . (5) 



3/2 



1/2 
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The decay width for the pseudoscalar A boson can be represented by exchanging exponents 3/2 <-> 
1/2 and m^o <-> m^o in Eq.(5). When m t + m c < M s < 2m t , the dominant decay modes of h° and 
A are h°, A — > cc, 66, tc + tc, whereas when M s > 2m t , the final state tt decay channel is open, and 
their decay widths are rather large due to the large masses of M s and m t . 

We denote 9 as the scattering angle between one of the gluons and the final top quark. Then we 
express all the four-momenta of the initial and final particles in the center-of-mass(CMS) by means 
of the total energy yf§ and the scattering angle 9. The four-momenta of top quark and charm quark 
are pi and p 2 respectively and are read 



(6) 



Pi = [e u ^Ej -m 2 sin9, 0, yj E% - rrftcosO^ , 
p 2 = (e c , —\J El — m 2 sin9, 0, —^E 2 — m 2 cos9) , 

where 

Et = \ (VI + K - ml) /VI) , E c = ^(V§-(m*- m^/VI) . (7) 
p 3 and p 4 are the four-momenta of the initial gluons and are expressed as 

Ps = Qv 7 !, o, o, ^Vtj , p 4 = (^v 7 !, o, o, . (8) 

The corresponding matrix element for all the diagrams in figure 1(a) and figure 1(b) is written as 

M = Tr(T a T b )5jiM Sl + (f abc T^)M^ + (T^ nl )M i + (T> m T^)M* (9) 

The upper indexes s±, §2, t and u represent the amplitudes corresponding to the s-channel dia- 
grams Fig. 1(a) (25~ 28), s-channel diagrams Fig. 1(a) (29 ~ 34) and Fig. 1(b) (13 ~ 15), t-channel 
and u-channel diagrams in figure 1(a) and figure 1(b) respectively. The T a (a =1 — 8) are the 377(3) c 
generators introduced by Gell-Mann and f a b c are the antisymmetric SU (3) c structure constants. The 
subscripts = 1,2,3) of T a represent the color of final state top quark and charm quark re- 

spectively. The variables s, t and u are usual Mandelstam variables in the center of mass system of 
gluon-gluon. Their definitions are: 

8 = (Pi +P2) 2 = (P3 +Pi) 2 , t = (pi -p 3 ) 2 = (p 2 -Pi) 2 , / 1q n 
U = (pi -Pi) 2 = (P2 -PS?- 

We collect all the explicit expressions of the amplitudes appearing in equation (9) in the appendix. 
The total cross section for gg — > tc + tc can be written in the form 

where \M\ 2 is the initial spin-averaged matrix element squared and = l/2(m 2 + m 2 — s) ± 
JE 2 — m 2 \f§. The cross section for pp — > gg — > tc + tc + X is conveniently written in terms of the 
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rapidities y\ and y2 of the two jets (finial states) and their common transverse momentum p?. Here 
we neglect the intrinsic transverse momentum carried by partons. It is 

j f°" j — = —r^-fg(xu Q 2 )f g (x2, Q 2 )°(gg -> tc + ic at 8 = rs), (12) 
dyidy2dpT s 

where y/s and \/I denote the proton-proton and gluon-gluon cm. energies respectively and s = sr. 
f g {xi, Q 2 ) is the distribution function of gluon in proton. 
Defining 

y* = ~(yi-y 2 ) (13) 
yboost = -ziyi + yi). (14) 



and 



We may write 



Ap 2 

T 

i Mi 



r coshV (15) 



xi = ^Tc Vhoos \ x 2 = \/re VboOBt . (16) 

In our numerical calculation we adopt the MRS set G parton distribution function f g (xi, Q 2 ) |T^| 
and let the factorization scale Q 2 = s. The numerical calculation is carried out around the Tevatron 
and LHC energy ranges. 

III. Numerical Results and Discussions 

In the numerical calculation we take the input parameters^ as = 4.5GeV, m c = 1.35GeV, 
m t = 175GeV, M w = 80.2226GeV A , G F = 1.166392 x 1(T 5 (GW)- 2 and a = 1/128. We adopt a 
simple one-loop formula for the running strong coupling constant a s as 

^(m) = - ~ 33- 2 n / : — r, f ^ n - ( 14 ) 
1 H Tr-t-otsVrnz) In ' — 



6-7T s v " I \ mz 

where a s (mz) = 0.117 and rif is the number of active flavors at energy scale /i. 

Figure 2 shows the cross sections for gg —>■ tc + ic as a function of the masses of the Higgs bosons 
M s . The cross sections are displayed with the three values of the gluon-gluon CMS energy 200 GeV, 
400 GeV and 500 GeV respectively. Because there is no stringent bound on the Higgs bosons masses, 
we choose M s in the range from 50 GeV to 800 GeV. The peak of each curve comes from s-channel 
resonance effects, where M s = m h o = m^o ~ y/B. From these curves we find that the cross section 
can be obviously enhanced when M s gets close to y/§. 

Figure 3 shows the cross sections of gg —>■ tc + ic as a function of and the three curves 
correspond to the M s values 100 GeV, 250 GeV and 500 GeV, respectively. For M s = lOOGeV, 
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the effects of the widths of the Higgs bosons are not obvious, and s-channel resonance effects are 
suppressed, since V§ is far beyond the Higgs boson masses M s . Therefore the curve of its cross 
section is relative flat with the increasing of y/S. When approaches the value of M s , such as 
M s = 250GeV and M s = 500GeV, the cross sections will be enhanced by the s-channel resonance 
effects, and the width effects become stronger, since the h°, A — > tc + tc channels are opened. In 
Fig. 3, we can see that the curve for M s = 500GeV shows a sharp peak around the position at 
\/l ~ 500GeV due to the s-channel resonance effects and large width effects of h° and A . 

In figure 4 and figure 5 we show the transverse momentum spectrum of the top quark at the 
Tevatron and LHC energies respectively, where we assume the off-line analysis will require at least 
the event selection criterion of involving one isolated high pr track with the cut of pseudorapidity 
\rj\ < 2. Again, the peaks on the curves of M s = 250 and M s = 500 show the s-channel resonance 
effects, where M s = m^o = m^o ~ \/§. 

In figure 6 and 7 we show the cross section of pp — > gg — > tc + tc + X as a function of center- 
of-mass energy of electron-positron system y/s. The cross section may reach 0.83 femtobarn when 
M s = 100 GeV and y/s = 2 TeV at the Tevatron and 131 femtobarn when M s = 100 GeV and 
y/s = 14 TeV at the LHC. For the Tevatron at 2 TeV we can expect about 4 raw events when 
M s = 100 GeV if we assume 5/6™ 1 integrated luminosity, and for the LHC at 14 TeV we can expect 
about 1.3 x 10 4 raw events if we assume 100/6 -1 integrated luminosity. Since the cross section 
of this process is roughly scaled by A 4 , if we let A ~ 1, the cross section will be 4 times larger. 
There are several potentially severe backgrounds from the SM to the signal. A top quark with a 
mass about 174GeV, decays dominantly to t — > W + b. In the tc — > W + bc — > l + vbc detection mode, 
the backgrounds are mainly from tt — > W + jets and tt WWbb — > l + vqq bb processes. The 
calculation shows that the cross section of top pair production will reach about 5 pb at the Tevatron 
for y/s = 1.8TeV and about 10 2 pb at the LHC for y/s = lATeV. If we use the tc + ic — > lepton+jets 
detection mode, the production rate of tc — > W + bc — > l + vbc can reach about 0.1 fb at the Tevatron 
and 14 fb at the LHC, while the possible background from tt — > WWbb — > l + vqq bb would be 
about 0.8 pb at the Tevatron and some dozens picobarn at the LHC. The reduction of the these 
backgrounds is possible through various kinematics cuts on the transverse energy, on the rapidity 
of jets and leptons, or involving b-tagging. Due to the very small production rate for the signal of 
pp — > tc + tc, it is not so easy to suppress its backgrounds. Therefore the further precise analyses are 
necessary to exclude these backgrounds. 

In summary, from our calculation, we can conclude that if we could efficiently exclude the severe 
backgrounds of the tcitc) production signal, it would be possible at the Tevatron and the LHC that 
the process pp — > gg — > tc + tc could be used to probe the flavor changing interactions in the context 
of the THDM III. ' 
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Appendix 

We adopt the same definitions of one-loop A, B, C and D integral functions as in Ref.|16j and 
the references therein. The dimension D = 4 — e. The integral functions are defined as 

A (m) = -M^ / d D q— -, 

{Bi, B^; B flu }(p,m 1 ,m 2 ) = — d q 

ITT J 



[q 2 — mf] [(q + p) 2 — m 2 



{C ; C p ; C^] C liup }(pi,p 2 ,m 1 ,m 2 ,m 3 ) = -^—^ 



21 ' 

— 7/1 

4-D 



x I d D q {1? 9ju> 9am/> ?A«/p} 



[q 2 - ml] [(q + px) 2 - m 2 .} [(q + p 1 + p 2 ) 2 - m§] ' 

(27T/i) 4_D 

{D Q ;D^D^;D^ p ;D^ pa }(p 1 ,p2,p 3 ,m 1 ,m 2 ,m 3 ,m i ) = — — 

x J d <?{1; q^t', q^iu'i Qfivp] Qfivpa} 

x {[q 2 - m\\ [(q + pif - m 2 2 } [(q + p 1 + p 2 f - 777,3] [(q + p± + p 2 + p 3 f - ml}}' 1 . 

In our calculation we take the strange quark mass m s = 0. The M Sl and M S2 in equation (9) can 
be written as 



M§1 = 16 ° 2 Af& mt V m t m c e » (P3) ^ (P4)U(P1 ) 

■{2a h om 2 (C + 4C 22 - 4C 23 )\p 3 , -p x - p 2 , m t , m u m^iplPx + P1P2 + P1P2 + P2P2) 
+2a h om 2 g^ u (B [-p 1 - p 2 , m t , m t ] - {{p x + p 2 ) ■ p 3 C + 4:C 24 )[p 3 , -p x - p 2 , m t , m t , m t ] 
+2ia A om 2 C \p 3 , -pi - p 2 , m t , m t , m t \e^ ua(3 'j 5 (p^p§ + p 2 p 3 )}v(p 2 ), 

(Al) 

where 

1 1 

a h° — ~ 2 — —rp, > a A° — ~ 2 — • 

s - mf i0 + xi hom h o s - m A o + u Aom A o 
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M ' S2 = l2Z^Mls m W m i m c e ^) e APi)u(pi) 

•{-16m t (C n - C12 + C 21 - C 23 ) + P2, M a , m t , m t ](p^ - pfr#) 

+ (8(2C 24 + m 2 t (C - C n + C 12 - C 21 - C 22 + 2C 23 ) - 2( Pl • p 2 )(C 22 + C 23 )) 

+p 2 ,M s ,m t ,m t ] +4m t 2 J Bi[-pi,m t ,M s ]- 
Am t m c Bi [p 2 , m t , M s \ + 4m^ Si [-pi, m 6 , M s ] + 
4m t m c Bi \p 2 , m b , M s \ + m^Bi [-p l ,m t , M s ] 
-m t m c B x \p 2 , m t , M s \ - m 2 B x \-p\,m b , M s ] 
+m c m t Bi[p 2 , m b , M S ]){^ V P\ ~ 7 M Pi + 1 V P2 ~ 1^P2 + 9^ fa) 
+ {-A(2m t C 2A - 2m 3 t C + m t {m 2 + 2( Pl ■ p 2 ) - 4(pi • p 3 ))(Cn - Ci 2 + C21 + C 22 
-2C 23 ) + A(p 2 ■ p 3 )(m c C 12 - m t C 22 + m t C 23 ))[- Pl , Pl + p 2 , M s , m t , m t ] 
-(m c m 2 C + 2m t C 24 + m*(C n - C 12 + C 21 + C 22 - 2C 23 ) + 2m t (( Pl ■ p 2 )- 
2(pi • p 2 ))(Cn - Ci 2 + C 2 i + C 22 - 2C 23 ) + 4(p 2 • p 3 )(m c C 12 
+m t (C 23 - C 22 ))[- Pl , Pl +p 2 ,M s ,m b ,m b ] - m%B 1 [-p 1 ,m t ,M 8 ]+ 
m 2 t m c B x \p 2 ,m u M^ - m^B 1 [-p l ,m b , M s ] + (m c m 2 + 
2m c p x ■ p 2 - Arricpx ■ p 3 )B 1 \p 2 , m b , M s ])g^ u - 4m t (C u 

-C12 + C 2 i - C 23 ) + P2, Af a , m 6 , m 6 ] (rfp 2 - PipZ + l 5 PiP 2 ~ I5P1P2) 

+ (2(2C 24 + m c m t C - m 2 {C lx - C 12 + C 21 + C 22 - 2C 23 ) - 2( Pl ■ p 2 )(C 22 - C 23 )) 

+p 2 ,M s ,m b ,m b } + 4m 2 B 1 [-p 1 ,m b , M a \- 
km c m t B x \p 2 , m b , M s })(YPi ~ YPi + YP2 ~ YP2 + 9*" fa) ~ (™ c m 2 C - 2m t C 24 
-m t 3 (C n - C12 + C 21 + C 22 - 2C 23 ) - 2m t ((p! • p 2 ) - 2( Pl ■ p 2 ))(C n - C 12 + C 21 
+C22 -2C 23 ) + 4(p 2 -p 3 )(m c C 12 -m t (C 23 - C 22 ))[-pi,Pi + P2, M s: m b: m b ]-f 5 g^ 
+ ((4C 24 - 2m c m t C - 2m\ (C n - C 12 + C 21 + C 22 - 2C 23 ) + 4(pi • p 2 )(C 23 
-C' 22 ))[-pi,pi +p 2 , M s ,m 6 ,m b ]) +2m^i[-pi,m 6 ,M s ]- 
2m c m t B l [p 2 , m b , M s ])>y 5 (YPi - YPi + Yp 2 - 1^P 2 - g^ u fa)}v(p 2 ), 

(A.2) 

The amplitude of M* can be written as 

■ 2 2 

M< = 64^M^ m * V m * m ^ M (P3 ) ^ (P4 ) U (Pi ) ( flPiPi + f2PlP2 + hPlP2 + UP2P2 

+hYp1 + hi^Pl + hi v rt + !^p v 2 + hi 11 !" + hoYY + hxfapIPl 

+h2fapip v 2 + fufartrt + fufapfrz + fi5farrf + hsfarpi + fnfarrt 
+fisfarP2 + fwfarY + f^fa-fr + Ai^pi + f^rta + r^iM <a 3) 

+/475^ + fhrfrf + r^P\ + fil^Pt + /s757 M P2 + /9757 M 7 iy 

+ /( 757 !y 7 M + fil-ysfarftt + f^faPlPl + fi 3 ^faPlP2 + fult>faP2P2 

+f[ 5 ^fal u P^ + f^faYPi + fulsfaYP* + flslsfa-fft 

+f[^fa^i u + f^farr)^), 

where the fiS and //s are expressed explicitly as, 

f[ = 4(m c (D 38 - D 310 ) +m t (D 32 - D 36 ))\p 3 , -p u -p 2 ,m b ,m b , M s ,m b ] (AA) 
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f[(m t -> -m t ) + 8m t (D u - D 12 + D 21 - D 24 - D 25 + D 2Q )\p u ~P3, ~P4, M s ,m t ,m t ,m t ] 

(A5) 



Aai{m c C 22 - m t C 12 - m t C 23 )[-p 4 ,p 2 ,m b ,m b , M s ] 

+4(m c (£> 39 - D310) - m t (D 22 - D 2& + £> 36 - D 38 ))\p 3 , -p u -p 2 , m b , m b , M s , m b \ 



(A6) 



/ 2 = f 2 (m t _> - mt ) + 8aim t Cn[-p 2 ,P4,M s ,m t ,m t ] ,^ _ 
-8m t (Di2 + -D 2 4--D26)[Pi,-P3,-P4,M s ,m t ,m t ,m t ] 

= 4(m c (-D 39 - L> 37 ) - m t (D 310 - D 38 ))[p 3 , -p u -p 2 ,m b ,m b , M s ,m b ] (A.8) 

h = f'^t -> -mt) ~Sm t (D 25 - D 2& )\p u -p 3 ,-p A ,M s ,m u m u m t ] (A 9) 

4(m c (D 33 - £> 37 ) + m t (-D 2 3 - L> 2 6 - Ano + Ak>))[P3, -p 2 ,m b ,m b , M s ,m b ] (AlO) 
f\ = f\(jn t -> -m t ) + 8m t D 2(i [p 1 ,-p 3 ,-p 4 ,M s ,mt,m t ,m t ] 



-2a 2 a 3 m%Bi[-p 1 ,m b , M s ] - 2a 1 a 2 (m 2 - 2{p x -p 3 ))-Bi[-Pi +p 3 ,m 6 , M s ] 

-2aia 3 m c m t Si[p 2 ,m 6 ,M s ] + 2a 2 (m c m t (C + 2C U + C 2 i) - 2C 24 )[pi, -p 3 , M s , ra t , m t ] 

+2a x (-2C 24 - m^Co + m2(C 22 - 2C 23 ) - m c m t C 12 + m?(Cn + C 2 i) 

+2(pi -p 2 -p 2 - p 3 )(C u - C12 + C 2 i - C 23 ) 

-2(pi •p 3 )(Cn + C 2 i))[-p 4 ,p 2 ,m 6 ,m 6 , M s ] 

+2(2^27 + 2/^3X2 + m 2 b D - m 2 c D 23 - m c m t (D 22 - D 26 ) 

+2(p 2 • p 3 )(D 25 - D 26 ))[p 3 , -pi, -p 2 ,m b ,m b ,M s ,m b ] 

(A12) 

/5 = /s^t ~ m <) + 4a 2 a 3 m^ [-pi, m u M s \ 

+Aaia 2 m 2 t B Q [-p 1 + p 3 , m t , M s \ - Aaia 3 m 2 t B [p 2 , m t , M s \ ,^ ^ 

-4a 2 m 2 (C + C 11 )\p 1 ,-p 3 ,M s ,m t ,m t \ - 4a 1 m 2 [-p 2 ,p 4 , M s ,m t ,m t ] 
+4m 2 (D + D n - £>i3)[pi, -P3, -p 4 ,M s ,m t ,m t ,m t ] 

= 2(4(^311 - ^312) + m^u - #12) - m*(£>37 - ^39) 

+m 2 (D 32 - D 36 ) + 2(pi • p 2 )(£>38 - Auo) - 2(pi • p 3 )(£>22 - D 24 - £34 + D m ) (A14) 
+2(pi -P2)(D 35 - D 310 ))[p 3 , -p 1: -p 2 ,m b ,m b ,M s ,m b } 

h = n (a 15) 



10 



= 2(-4D 3 i3 - m 2 b D 13 + m 2 c D 33 + m c m t D 23 + m 2 (D 26 + D 38 ) + 2( Pl ■ p 2 )(D 23 + D 39 ) 
-2(pi -P3)(D 25 + D 310 + D 23 + D 37 ))\p 3 , -pi, -p 2 ,m b ,m b ,M s ,m b ] 

(A16) 

/ 7 = /^(mt -> -m t ) - Am 2 t Di 3 [pi, -p 3 , -p 4 , M s , m u m t , m t ] (A17) 

f% = -4:a 1 ( Pl - P3 )(C 1 2 + C 23 )[- P 4,P2,m b ,m b ,M s \ 

+2(2(D 27 + 2D 311 + D 313 ) + m 2 (D + D n ) - m 2 c (D 23 + D 37 ) 

-m 2 (D 22 + D 36 )\p 3 , -p u -p 2 ,m b ,m b ,M s ,m b ] (A18) 

-(Pi ' P 2 )(^26 + Auo) + 2(pi • p 3 )( J Dl2 - D 13 + 4D 24 - ^26 + D 34 ) 

+2(p 2 -p 3 )(D 25 + D 35 ))\p 3 , -pi, -p 2 ,m b ,m b ,M s ,m b ] 

h = fs (A 19) 



= -a 2 a 3 mfS 1 [-pi,m 6 ,M s ] 

+2aia 2 m c (pi -p 3 ).Bi[-pi + P3,m b , M s ] - a 2 m t (2C 2A 

-m 2 {C + 2C U + C 21 ) + 2(p! • p 3 ){C + C n + C 12 + C 23 ))[p!, -p 3 , M s , m t , m t ] 
-2aim c (pi •p^Ci^-^p^m^m^M,;] + (4m c (0 3 i 3 
+m2 J D 13 - m^s) + m t (2D 27 + 4D 312 + mg(L» + D 12 ) - m 2 c (D 23 + D 39 ) 
-m t (m c D 3S - m t D 22 - m t D 32 )) - 2(p x ■ p 2 )(m c D 39 

+m t D 26 + m t D 38 ) + 2{p ± ■ p 3 )(m c D 12 - m c D 13 + m c D 310 + m t D 24 + m t D m ) 
+2(^2 -P3)(m c D 37 + m t D 25 + m t D 3W ))[p 3 , -p u -p 2 ,m b ,m b , M s ,m b ] 

h = f'%{ m t -»■ -m*) - 2a 2 a 3 m t 3 S [-pi, m t , M s ] 
-4aia 2 m t (pi -p 3 )S [-pi +p 3 ,mi,M s ] 
+Aa 2 m t (pi ■ p 3 )Cobi, ~P3, M s ,m t , m t ] 
+Aaiin t (pi ■ P3)C [-p 2 ,p 4 , M s ,m t ,m t ] 
+4m t (D 27 + D 3 n - D 3 i 3 )[pi, -p 3 , -p 4 , M s ,m t ,m t ,m t ] 

f[ = -2(m c D 313 + m t .D 3 i 2 )[p 3 , -p 2 , m 6 , ra 6 , M s , ra 6 ] 



/io = f[ ( m t -> -"it) - 4m t ^ 27 [pi, -p 3 , -p 4 , M s , mt, m t , m t ] (A23) 

/n = 4(L> 22 - L> 24 - L> 34 + D 36 )[p 3 , -p u -p 2 , m b , m b , M s , m b ] (A24) 

/n = fn (A25) 

/( 2 = 4(L> 13 - D l2 - 2D 24 + 2D 26 + D 3l0 - D u )\p 3 , - Pl , -p 2 , m b , m b , M s , m b ] (A26) 
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(A20) 

(A21) 
(A22) 



fn = f[ 2 (A.27) 

fn (A28) 

fa = fiz (A29) 

/( 4 = 4(^23 - ^25 - D 35 + £>37)[P3, ~P2, m b , m b , M s , m b ] (AM) 

fu = fu (A31) 

f[ 5 = -2a 2 m t (Ci 1 -Ci2 + C 2 i- C 2 3)[pi,-P3,M s ,m t ,m t ] (A 32) 
+2(m c (D 13 - D 12 + D 26 ) + m t (D 22 - D 24 ))\p 3 , -p x , -p 2 ,m b ,m b , M s ,m b ] 

fib = f[ 5 ( m t -> ~ m t) + 4a 2 mtCii[pi, -p 3 , M s , m t , m t ] ^ 
-Am t (D u - -Di 3 )[pi, -p 3 , -pA,M s ,m u m t ,m t \ 

f[ 6 = 2(m c (D 25 - D 26 ) -m t (D 22 - D 24 ))[p 3 , -p u -p 2 ,m b ,m b , M s ,m b ] (AM) 

fia = fia( m t -»■ ~ m t) +4m t (D n - D 12 )\p u -p 3 , -p A ,M s ,m t ,m u m t \ (A35) 

f 17 = 2(m c D 23 - mt(D 25 - D 26 ))\p 3 , -p u -p 2 ,m b ,m b , M s ,m b ] (A.36) 

fn = fn(mt -> -m t ) + ^m t D 13 \p u -p 3 , -p 4 , M s , m t , m t , m t ] (A.37) 



f[ 8 = 2a l (m c C 22 -m t Ci 2 -m t C 23 )[-p i ,p 2 ,m b ,m b) M s ] M 38) 

+2(m c (D 2 5 - A23) + m t (D 12 - D 13 + D 24 - D 26 )) [p 3 , -p u -p 2 , m b , m b , M s , m b \ 

fis = /{ 8 ( m * - m t) +4:a 1 m t C 11 [-p 2 ,p4,M s ,m t ,m t \ (A 39) 

-4:m t D 12 \p 1: -p 3: -p4,M S: mt : mt : mt\ 
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f[g = -a 2 a 3 m 2 t B 1 [-p 1 ,m b ,M s } 

-a x a 2 (m 2 t - 2(pi -p 3 )).Bi[-pi + p 3 ,m b , M s \ 
-a x a 3 m c m t B x \p 2 , m b , M s \ - a 2 (2C 24 - m t (m c C + m t C u 
+m t C 21 ) + 2(Ci2 + C 23 )(pi -p 3 ))[pi, -p 3 ,M s ,m t ,m t ] 
+a 1 (-2C 24 - mgC + m2(C 22 - 2C 23 ) - m c m t C X2 
+m 2 {C u + C 2 i) + 2((pi • p 2 ) - (pa • p 3 ))(Cn - C 12 + C 21 - C 23 ) 
-2(pi -p 3 )(Cii + C 2 i))[-p 4 ,p 2 ,m b ,m b , M s ] 

+(4(D 27 + Am) + 2mg(A, + D n ) - 2m 2 c (D 23 + £37) - D 13 m c m t 
-m 2 (D 12 + D 22 + ^e) - 2( Pl -p 2 )(D 12 + 2D 26 + £> 310 ) + 2( Pl • p 3 )(D 12 
+2D 24 + D 34 ) + 2(p 2 ■ p 3 )(D 13 + 2D 25 + D 35 ))\p 3 , -p u -p 2 , m b , m b , M s , m b ] 



where 



1 1 A 1 

ai = j -, a 2 = j and a 3 = — = -. 

t — rrit t — m 2 mf - m z c 

M" = M* {p 3 ^p A , //<-»• v, i <-> u) (AAA) 



(AAO) 



f[ 9 (mt -> -m t ) + 2a 2 a 3 m 2 B [-pi, m t , M s \ 

+2a 1 a 2 m 2 B [-p 1 + p 3 , m t , M s \ - 2a 1 a 2 m 2 B \p 2 , m t , M s \ ,^ ^ 

-2a 2 m 2 C \pi, -p 3 , M s , m t , m t ] - 2a 1 m 2 C [-p 2 ,P4, M s , m t , m t \ 
+2m 2 D \pi, -p 3 , -pa, M s , m u m t , m t ] 

/20 = -2D 311 \p 3 , -pi, -p 2 ,m b ,m b , M s ,m b ] (AA2) 
/20 = f 20 (A43) 
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Figure Captions 

Fig.l(a)(b) The Feynman diagrams of the subprocess gg — > tc. 

Fig. 2 Total cross sections of the subprocess gg — > tc + tc as function of M s . The solid curve is 
for Vs = 200GeV, the dashed curve is for \f§ = AOOGeV and the dotted curve is for Vs = 500GeV. 

Fig. 3 Total cross sections of the subprocess gg — > tc + ic as function of V§. The solid curve is for 
M s = 100 GeV, the dashed curve is for M s = 250 GeV and the dotted curve is for M s = 500 GeV. 

Fig. 4 The transverse momentum spectrum of top quark for the Tevatron at 2 TeV. The solid 
curve is for M s = 100 GeV, the dashed curve is for M s = 250 GeV and the dotted curve is for 
M s = 500 GeV. 

Fig. 5 The transverse momentum spectrum of top quark for the LHC at 14 TeV. The solid curve is 
for M s = 100 GeV, the dashed curve is for M s = 250 GeV and the dotted curve is for M s = 500 GeV. 

Fig. 6 Total cross sections of the process pp — > gg — > tc + tc as function of y/s at the Tevatron. 
The solid curve is for M s = 100 GeV , the dashed curve is for M s = 250 GeV and the dotted curve 
is for M s = 500 GeV. 

Fig. 7 Total cross sections of the process pp — > gg — > tc + tc as function of y^s at the LHC. The 
solid curve is for M s = 100 GeV, the dashed curve is for M s = 250 GeV and the dotted curve is for 
M s = 500 GeV. 
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